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1 Introduction 



The investigation of the renormalization-group (RG) properties of composite oper- 
ators in QCD is an important ingredient in the numerous applications of QCD in 
particle physics. In the present paper we shall consider some aspects of the RG 
properties of heavy baryonic currents in the framework of the Heavy Quark Effec- 
tive Theory (HQET) [1], [|, in particular we consider the anomalous dimensions of 
the currents associated with the ground-state heavy baryons. It is well known that 
the heavy mass dependence of operators and matrix elements in QCD can be very 
well organized by employing an 1/mQ-expansion provided by HQET. Here we work 
to lowest order in the heavy mass expansion, i.e. we treat the heavy quark fields as 
effective static fields. While this approximation should be a good starting point to 
understand the properties of heavy hadrons containing one heavy quark (b or c) it 
leads to significant simplifications in the calculation of the Wilson coefficients that 
appear in the Operator Product Expansion (OPE) and in the calculation of RG 
coefficients such as the anomalous dimension coefficient dealt with in this paper. 

The anomalous dimensions of composite operators are basic ingredients in all 
applications of the OPE method in QCD such as QCD sum rules || and related 
approaches. In this paper we undertake to calculate the two-loop anomalous di- 
mensions of heavy baryon currents in the static approximation. The anomalous 
dimensions of the baryonic currents are one of the ingredients needed in the analysis 
of the correlators of two baryonic currents via their short distance expansion. The 
ensuing QCD sum rules allow one to compute the mass of heavy baryons and their 
residues in terms of basic non-perturbative QCD parameters. 

Another important application where the anomalous dimensions enter is the de- 
termination of the universal non-perturbative baryonic Isgur-Wise function £(u;) [§j 
from three-point QCD sum rules. Of particular importance is the calculation of the 
slope of the Isgur-Wise function at zero recoil. Knowledge of the two-loop anoma- 
lous dimensions is important when one discusses the matching of HQET baryonic 
currents with the corresponding currents in full QCD. Furthermore, the anomalous 
dimensions are needed for the calculation of some basic QCD matrix elements that 
appear in QCD sum rules, as for example the expectation value of the kinetic energy 
//^, the expectation value of the chromomagnetic energy /jLq, the spin-dependent axial 
current matrix element the couplings of heavy baryons with pions and photons, 
and the magnetic and axial moments of heavy baryons. 

The general structure of heavy baryon currents has the form (see e.g. [|[ and 
refs. therein) 

J = [ q * T CTTq3]T'Q k e ljk . (1) 

Here the index T means transposition, C is the charge conjugation matrix with the 
properties C^JC -1 = — 7 M and C^C~ l = 75, i,j,k are colour indices and r is a 
matrix in flavour space. The effective static field of the heavy quark is denoted by Q. 
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The currents of the Ag and the heavy quark spin doublet {Sq, Sq} are associated 
with the spin-parity quantum numbers j p = + and j p = 1 + for the light diquark 
system, respectively. For each of the ground state baryon currents there are two 
independent current components J\ and J r 2 . They are given by |5j, [| 

J ax = [q iT CT-f 5 q j }Q k e ijk , J A2 = [q iT CT^ q j }Q k e ijk , 
Jvi = [q iT Cr^q j ] ■ 775<2%'fc, ^S2 = [q lT Cr^q j ] ■ 775<2%'fc, 
4i = [q iT Cr^}Q% k + ^[q lT Cr^}^Q k e ijk , (2) 

where J^n and Js. 2 satisfy the spin-3/2 condition 7 Js*i = (i = 1, 2). The flavour 
matrix r is antisymmetric for Aq and symmetric for the heavy quark spin doublet 
{Sq,Sq}. The currents written down in Eq. (^]) are rest frame currents. The 
corresponding expressions in a general frame moving with velocity four-vector v fl 
can be obtained by the substitutions 7o — > t and 7 — > 7^ = 7 M — M . 

The one- loop renormalization of these currents was considered in . The renor- 
malization of the corresponding light baryonic currents to one-loop order was pre- 
sented in [£| and to two- loop order in 0. In this paper we determine the two-loop 
anomalous dimensions of the heavy baryon currents. We evaluate the numerous 
two-loop diagrams with the help of REDUCE and MATHEMATICA packages writ- 
ten by us, which are based on the recurrence-relations in QCD and HQET presented 
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The paper is organized as follows. In Sec. 2 we furnish some general remarks 
on anomalous dimensions, discuss the one-loop renormalization of the general heavy 
baryon currents ([!]) and give results on the one-loop anomalous dimensions in a 
general covariant gauge. In Sec. 3 we discuss general features of the two-loop renor- 
malization procedure and present our general two-loop results, where we specify to 
the Feynman gauge. In Sec. 3 we also derive a powerful consistency check on our 
calculation, and we compare our results with previously published two-loop results 
in the mesonic sector. In Sec. 4 we consider applications of our general results to 
the ground-state heavy baryon currents (0). We work in the MS-renormalization 
scheme throughout. As concerns the treatment of 75 we give our results for two 
different schemes, namely the 't Hooft-Veltman 75-scheme and the naive anticom- 
muting 75-scheme. Finally, we show how the results in the two 75-schemes become 
related through a finite renormalization prescription as has been demonstrated pre- 



viously in a different context by Trueman and Larin [13|, |14|| . Sec. 5 contains our 
conclusions. 



3 



2 Renormalization procedure and one loop result 



The renormalization procedure consists in the renormalization of all bare UV-diver- 
gent operators On in terms of a renormalization factor Z according to 



On = ZO. 



(3) 



Using dimensional regularization in D = 4 — 2e dimensions and the MS-scheme Jl5| , 
the factor Z is a polynomial in inverse powers of e. The Z-factor is constructed in 
such a way that the operator O is finite in physical four- dimensional space. After 
renormalization, the operator O can be seen to depend on the subtraction scale 
(normalization point) fi. The anomalous dimension of the operator O is then defined 
by 

din Z(a([i), a; e) 
7 = rflm^O ' 
where a is the renormalized gauge parameter in the general covariant gauge (with 
a gluon propagator proportional to —g^ v + (1 — a)k^k v /k 2 ) and a(/i) is the renor- 
malized coupling constant in four- dimensional space. The Z-factor and thereby the 
anomalous dimension 7 is of course specific to the operator which it renormalizes. 
Thus they would carry an additional index which we suppress for the moment. The 
unrenormalized gauge parameter a and the unrenormalized coupling constant an 
defined in D dimensions are related to the corresponding renormalized quantities a 
and a(/j,) in four dimensions by 



«o = a(/x)/i e Z a {a{jj), a; e), clq = aZs(a(n) } a; e). 



One-loop MS-results for the factors Z a and Z 3 have been given e.g. in |L6 
read 



Z a 
Z, 



1 + 



a s 

47T6 

4tc€ 



—C a - -T F N F 
3 3 
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6 3 . 



(5) 
They 

(6) 
(7) 



Note that there is a sign difference in the second 0(a s ) terms in Eqs. (|6|) and (0) 



relative to that in |16[ due to a difference in the sign convention for e. The factors Z a 
and Z3 determine two universal RG functions (3 and 5 which denote the anomalous 
dimensions of the coupling and gauge parameters, i.e. 



/3(a(/i), a) 



dlia Z a (a(fi), a; e) 
dln(fi) 



and £(a(/i),a) 



dhi Z 3 (a(/i), a; e) 
dlia(fi) 



(8) 



In perturbation theory the renormalization factor or function Z(a(/i), a; e) will be a 
double series in the renormalized coupling constant a s = a(/x) and in inverse powers 
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of e. We therefore write the double series expansion 

m=l fc=l v ^ ;l 7 fc fc=l fc 

In contradistinction to this the anomalous dimension 7, as well as the functions (3 
and 5, are only expanded in powers of the coupling constant a s since they are finite 
quantities for D — > 4. One has 
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E fS 7 (m) , /?=E& *=Efe do) 



m=l ^ m^l V47F ^ m^l V4?r 



Using these functions and differentiating Z(a(/i), a; e) with respect to ln(/x) one can 
compare coefficients at different orders of l/e k ||16|| . For /c = one obtains 



7 = -2^- rn) 

a In a s 

and for k > 

-2^±L=f 7 -/J Jn »— i3 A-)z 4 . (12) 
a m a s y In a s alna/ 

The relations (0) and fll2|) allow one to determine the anomalous dimension on the 
one hand and, on the other hand, provide for a connection between the coefficients 
Z m ^ at different powers of 1/e which allows one to check on the consistency of the 
two-loop calculations. For example, the first order contributions in 1/e determine 
the one- and two-loop anomalous dimensions 

7 {1) = -2Z hl and 7 (2) = -4Z 2)1 . (13) 

For the two-loop coefficients, Eq. (O) leads to the consistency condition 

- 4Z 2 , 2 = (V } - ~ ^g^) (14) 
2.1 One- loop results 

Let us first consider the one-loop renormalization of the effective heavy baryon 
currents ([]]). The bare light quark and bare effective heavy quark fields are related 
to the renormalized fields by 



ZWq, Q = zfQ. (15) 



In the MS-scheme with D = 4 — 2e space-dimensions we have 
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where go is the bare coupling in QCD and do is the bare gauge parameter. We use the 
usual definitions for SU(N), i.e. C F = {N 2 C - 1)/2N C , C A = N c , C B = (N c + l)/2N c , 
and Tp = 1/2 for N c = 3, Np being the number of light quarks. The bare current 
is renormalized by the factor Zj, i.e. 

J = ^CTrq )T'Q = Z q Z% 2 Z v J = ZjJ. (17) 

As can be seen from Eq. fll7|) , the anomalous dimension of the full current J is 
composed of three parts associated with the renormalization of the light and heavy 
quark fields, and the renormalization of the vertex. One thus has 

7j = 27 9 + 1q + 7v- (18) 

The one-loop anomalous dimensions 7^ and 7q can be read off Eqs. flT5|) and ( [T5D 
and are 

^=C F a, 7 « = CV(a-3). (19) 
A concerns the corrections to the vertex (g r CTrgo)r'Qo, one obtains 

Z v = 1 + ^-{{n - 2) 2 + 3a - 1). (20) 
Aire 

The one-loop vertex Z-factor has been calculated in the general covariant gauge 
a / 1 since the consistency condition in Eq. ([14]) requires knowledge of the depen- 
dence of Zi t i on the gauge parameter a. In deriving Eq. fl2"0p we have parameterized 
the general light-side Dirac structure T in terms of antisymmetrized products of n 
Dirac matrices, Y = 7^ • • -7^"!. For the time being we assume that there is no 75 
in the Dirac string T and postpone the discussion on how to define 75 in D 7^ 4 di- 
mensions to the end of the paper when we present explicit values for the anomalous 
dimensions. The above form for the Dirac string T introduces a n- and s-dependence 
in the vertex corrections due to the identities 

la r la = hv = {-i) n {D - 2n)r, 7 r 7o = {-i) n s v. (21) 

Using the results of Eq. (|TB|) we immediately obtain the one-loop anomalous dimen- 
sions of the irreducible vertex 

t£> = -2C B ((n-2) 2 + 3a- 1). (22) 



According to Eqs. ([TJ) and (0), the one-loop anomalous dimension of the baryonic 
current is given by 

7J = T ( ~ 2 ° B ^ n ~ 2 ) 2 + 3a - 1) + 3C V(« - 1)) + 0{a 2 s ). (23) 

Note that the dependence on the gauge parameter a drops out of Eq. (p3|) when 
substituting Cp = 4/3 and Cb = 2/3. This is a necessary requirement since the 
anomalous dimension is a gauge-independent quantity. When one specifies to the 
Feynman gauge (a = 1), the results in Eqs. (|20|) , J2"2|) and (p3|) agree with the 
corresponding results presented in B. 
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3 Two-loop renormalization 



In this section we consider the two-loop renormalization of the heavy baryon cur- 
rent (|I|) in the MS-scheme where we now work in the Feynman gauge to simplify 
the calculation. First of all, the two-loop anomalous dimensions of the quark fields 



are well known and can be taken e.g. from [11, 17, 18, 19, pOfl. They are given by 



7f = C F (^C A - 2T F N F - 3 -C F ) , 7< 2) = C F (-fc A + fT F N F ) . (24) 

In order to determine the vertex renormalization factor Zy that renormalizes the 
bare proper vertex according to V = ZyY we need to compute the bare proper 
vertex Vo i n one-loop and two-loop order. The Dirac structure of the unrenormalized 
and renormalized vertices Vo an d V is given by the Born-term Dirac structure of 
the baryonic currents, i.e. by T <g) V in the notation of Eq. (|T]) (all heavy baryon 
currents are renormalized multiplicatively in the effective theory!). We can thus 
write V = K)F <8> r" and V = VT <g> V, where V is finite. By requiring V = Vo/Zy 
to be finite one can calculate Zy in terms of Vq. For this purpose one first evaluates 
the one- and two-loop diagrams in terms of the bare coupling constant go and, in 
case of the one-loop diagrams, in terms of the bare gauge parameter a^. One then 
substitutes for these in terms of the renormalized g s (gl = 47ra s ) and a with the 
help of the Eqs. (|5|), (^) and (|7]). The resulting expression for the bare proper vertex 
Vo may again be presented as a double series in a s and e, 

^o=££fe ? V m , k . (25) 

m=0fc=0 7 fc 



Thus the coefficients V m ^ and Z m ^ become related. For example, for the first few 
coefficients one has the relations 

#1,1 = ^,1, ^2,2 = ^2,2, Z 2 ,i = V 2 ,l - Vl,lVl, . (26) 

Note that in the relation for the Z-factor Z 2) i one also has to include one-loop 
contributions. 

Three different contributions to the two-loop vertex corrections Vo may be iden- 
tified. First one has the set of two-loop graphs where the two-loop contributions 
are associated with one of the heavy-light subsystems. Then there is the subset of 
two-loop graphs that are associated with the light-light system. Finally there are the 
irreducible contributions where the two loops connect all three quark lines. Hence 
we write 

V = 2V {hl) + V {ll) + Vt\ (27) 

where the factor of two on the r.h.s. accounts for the fact that there are two possible 
heavy-light configurations. The three types of contributions will be considered in 
turn in the following subsections. 
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3.1 The heavy-light system 

In this subsection we consider composite operators (qQ) with one massless quark 
field q and one effective static heavy quark field Q. There are 11 different two-loop 
diagrams contributing to the composite operator (qQ) which are shown in Fig. 1. We 
choose the momentum of the light quark to be zero and regularize all IR-divergencies 
by taking the heavy quark off its mass-shell. The bare proper vertex of the 
heavy-light system is calculated to two-loop order in the Feynman gauge (a = 1) 
using the algorithm developed in ||11|| . The result including the Born-term and 
one-loop contributions can be presented as 

Vi hl) = 1 + hC^ l) + £ AjCifr? , (28) 

i=l 

where we use the notation A/ := (g^/ (Att) d ^ 2 )(— 2oj) D ~ i . The contributions of the 
various diagrams in Fig. 1 are separately identified according to i = for the one- 
loop contribution and i = 1, . . . , 11 for the two- loop contributions. The associated 
colour factors are denoted by Cj, 

Co — Cb, C\ = C B , C*2 = Cb(Cb — \Ca), 

C 3 = C 4 = CbCf, C 5 = Cq = Cb{Cf — \Ca), (29) 

CV = CsNpTp, C§ = Cii = CbCa, Cq = Cio = \CbCa- 

One should keep in mind that the colour factors have to be calculated for the 
completely antisymmetric baryonic colour configuration q l q^Q k £ijk- Because of the 
difference in colour structure it is not a priori clear how the baryonic heavy-light 
two-loop results can be related to the mesonic case with colour structure qiQ^5j. 
However, it is not difficult to see that the corresponding colour-singlet colour factors 
for the mesonic case can be obtained from Eq. ( p9|) by the substitution Cb — » Cp. 
Note, though, that this argument cannot be turned around. The contributions to 
the baryonic anomalous dimensions cannot be read off from the mesonic ones. 

The coefficients are listed in Appendix A. The coefficients V^k are then 
given by 

Vif = C B a, Vif = 0, (30) 
V 2 {h 2 l) = C b (\c b -Ca), V 2 { f = -C B (C B (1-4C(2))-C A (1-C(2))). 

Using the relations ( ^6l) one can then calculate the coefficients Z n ^ which determine 
the two-loop anomalous dimension of the subset of the heavy-light graphs (with 
antisymmetric colour structure), 

7$) = Ca(4 " 16C(2)) - C B C A (4 - 4C(2)). (31) 
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There are two checks on our calculation. First, the coefficients ^2,2 and Z\^_ can be 
seen to satisfy the consistency condition Eq. (14). Second, as remarked on above, 
it is not difficult to transcribe the calculation to the mesonic case (colour-singlet 
heavy-light quark-antiquark system (jiQ^Sj with colour factors given after Eq. (|29|) ). 
After transcription to the mesonic case we find agreement of our results with the 



results given in [11, 19 



3.2 The light-light system 

The most difficult part of the calculation is the evaluation of the two-loop bare proper 
vertex of the light-light system. There are again 11 different Feynman diagrams 
corresponding to the two-loop corrections of the vertex with two light quarks (see 
Fig. 1 with the heavy quark line substituted by a light quark line). 

Again we need to compute the bare proper vertex V at one- and two-loop order. 
The external momentum of one of the light quark legs is set to zero. Using the power 
counting method, one can convince oneself that this procedure does not introduce 
any unwanted infrared singularities. All infrared singularities are regularized by 
choosing a nonzero momentum for the second light quark leg. On the other hand, 
setting one external quark momentum to zero, the ultraviolet properties of the 
diagrams are not affected. The advantage is that all diagrams now become two- 
point functions, which can be treated by using the algorithm presented in 0. In 
the first stage the one- and two-loop integrals are evaluated independently of the 
specific Dirac structure of the Dirac string T. In the second stage one specifies T and 
the general result can then be written down using the identities (pTD . The results of 
the calculation are again presented in coefficient form according to 

Vi u) = 1 + E ( XaCohSS + E A^C') ( S ( D " 2n )y> ( 32 ) 

j=0 \ i=l / 

where Xg ■— {g^/ {Art) D l 2 )(— l/p 2 ) D / 2 ~ 2 . The colour factors Ci have been calculated 
in Eq. (^), the index i is defined as in Eq. (f2S|). 

We need to explain the index j and the factor (s(D — 2n)y appearing in Eq. (|52"|). 
In the calculation of the one- and two-loop contributions one encounters different 
contracted forms of the general light Dirac string structure T according to 

T = T, Ti = 7 M 7 r7o7 At , r 3 = J^^p^JolpluJ^, (33) 

T 2 = 7p7u^7„7p, T 4 = 'Jp'Julpla^'lalp'lulp- 

Again, we defer the discussion on 75 to Sec. 4 assuming T to contain no 75. From 
the identities in Eq. (|2l|) one obtains Tj = (s(D — 2n)yT leading to the factor 
(s(D-2n)y in Eq. fl33|). 
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The coefficients bfj are listed in Appendix B. Then, organizing the contributions 



according to the double expansion in Eq. ( j25[) one has 

V}f = C B ((n-2) 2 + a-l), (34) 

Vif = C B ((n-2) 2 + (n-2)(s + 2) + a-l), (35) 
C 

V$ = -f (9C B (13n 4 -96n 3 + 248n 2 - 256n + 88 + 4(n-2) 3 s) 



36 

-72CV(n - 3) (n - 1) - C A (18n 2 - 144n 3 + 289n 2 + 128n - 444) 
-AN F T F (n 2 - 16n + 24)), (36) 



VS = — (3C B (n~2) 4 -C A (lln 2 -AAn + 39)+ AN F T F (n-3)(n-l)) . (37) 
g v / 

As in Sec. 3.1, the one-loop results have been calculated in the general covariant 
gauge whereas the two-loop results are given only for the Feynman gauge. The 
result for the two-loop anomalous dimension of the proper light-light vertex is then 
given by 

C 

7 gj = -jp ( - 9C B (5n 4 - 40n 3 + 104n 2 - 96n + 24 + 4(n - 2) 3 s) 

+72C F (n - 3)(n - 1) + C A (18n 2 - 144n 3 + 289n 2 + 128n - 444) 
+AN F T F (n 2 - 16n + 24)) . (38) 



We have checked that the coefficients and z[i satisfy the consistency condition 
Eq. ( |i~4| ) for any value of n. 

A further check of our calculation is provided by the colour singlet projection 
of the proper vertex for the light-light mesonic configuration (q l qjSj) which can be 
obtained from our result with the substitutions q T C — > q and C B — > C F . Collecting 
the Cir-contributions one has 

ft 

7 J 2 j = _Z(_9C Fn ( n _4)(5 ra 2 _20n + 16) 

+C A (18n 4 - UAn 3 + 289n 2 + 128n - 444) (39) 
+AN F T F {n 2 - 16n + 24)" 



After appending the anomalous dimensions 2^^ of the light quark fields one can 
then compare with previously published results in the mesonic sector. For n = 
Eq. (B^) coincides with the known anomalous dimension of the scalar current first 



written down by Tarrach |2~I |. For n — 1 one reproduces the well-known vanishing 
of the anomalous dimension of the vector current. For n = 3 one obtains agree- 
ment with the axial current case derived by Trueman ||13||, and, finally, for n = 4 
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with the two-loop anomalous dimensions of the pseudo-scalar current (with the 
't Hooft-Veltman convention for 75 to be discussed later on) which was calculated 
by Larin ||14|| . The anomalous dimension of the tensor current (with T = a^ u ; n = 2) 
was recently obtained by Broadhurst and Grozin J22J in an on-shell QCD calculation, 



who had also checked on the previously published results by considering a general 
Dirac structure for T. Again there is agreement. 

3.3 The heavy-light-light irreducible vertex 

Next one needs to calculate the three-quark irreducible proper vertex Vq . In this 
case there are altogether 8 diagrams. Four of them are shown in Fig. 2, the other 
four can be obtained by reflection on the heavy quark line which amounts to an 
exchange of the two light quark legs. 

To start with, each quark leg carries a momentum pi and thus every diagram 
depends on three external momenta. In diagrams (1), (2) and (4) of Fig. 2 we can 
safely set p\ = P2 = without introducing any infrared singularities, since they are 
naturally regularized by the off-shell energy u of the heavy quark. In diagram (3) 
we again set P2 = 0. Power counting arguments imply that the contribution of dia- 
gram (3) contains a term proportional to \n(p 2 ). We thus have to keep p\ finite. The 
integrals in diagram (3) cannot entirely be reduced to two-point two-loop integrals, 
as is the case for diagrams (1), (2) and (4). In fact the general structure of the 
two-loop diagram (3) is given by 

r dkdi (i- Pl rir (k-ty 

3 ^ up J {2ir) 2D {I - Pl ) 2 l 2 k 2 {u + kv) (k - I) 2 ' 1 ' 

where T M ^ p represents a string of Dirac matrices. The simplest way to evaluate 
this integral is to add and subtract an integral which has the same structure as 
the integral shown in Eq. (fEJ) but with the simplified propagator k v jk 2 instead of 
(k — l) u j(k — I) 2 . The simplified integral can be easily calculated even for pi ^ 
because it factorizes into two one-loop type integrals. The advantage is that the 
difference of the two integrals with propagators (k — l) u /{k — I) 2 and k v jk 2 is IR- 
finite as p\ — > and thus one can set p\ = when evaluating the difference integral. 
The latter integral can then be solved by the standard methods. The results of our 
calculation are again presented in coefficient form, 

Vt ] = 2 E ( E A^y) (s(D - 2n)Y, (41) 

j=0 \i=l ) 

where the factor of 2 accounts for the contributions of the remaining four reflected 
diagrams not included in Fig. 2. For the definition of the index j we refer to Eqs. (|32"D 
and fl3"3]). The index i identifies the relevant Feynman diagram in Fig. 2. The 
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coefficients are listed in the Appendix C. For the Z-factors and for the anomalous 
dimension we obtain 

%2% = C 2 B (2n 2 -8n + 9), (42) 
7( ( 2) = "4^i r) = -2CS(9n-10 + 6a)(n-2). 

It is interesting to note that the unrenormalized irreducible two-loop vertex contains 
a 1/ e 2 -contribution even though there is no corresponding one-loop contribution to 
the irreducible three-quark vertex. The coefficient of the 1 / e 2 -contribution is related 
to the 1/e-pole in the total one-loop result by the non-linearity in the interplay 
between Eq. (|i~3| ) and Eq. (|14|). We have added this comment because this interplay 
of one- and two-loop contributions is a feature specific to the baryonic case and does 
not occur in the mesonic case. The result in Eq. fl42| ) completes our calculations. 

4 Anomalous dimensions of Aq and Eg baryons 

We are now in the position to present our results for the anomalous dimension of 
the baryonic currents in Eqs. (@). We first sum up the results for the light-light, the 
heavy-light and the irreducible heavy-light-light vertex and obtain 

C 

7 { / 2) (n,s) = — (-9C B (5n 4 -40n 3 + 106n 2 -104n + 24 + 4(n-2)s + 32C(2)) 
9 ^ 

+C A (18n i - 144n 3 + 289n 2 + 128n - 516 + 72((2)) 

+72C F (n - 3)(n - 1) + AN F T F (n 2 - 16n + 24)). (43) 



To obtain the full two-loop result for the anomalous dimension one has to add the 
anomalous dimensions of the heavy and light quark fields from Eq. (|24D according 
to Eq. flT|). 

Up to this point we have steered clear of the 75-issue because we did not want to 
interrupt the general flow of arguments in the preceding sections. When one wants 
to apply the general result Eq. flilf ) to the Aq case in Eq. (H) one faces the notorious 
problem of how to generalize 75 to D 7^ 4 dimensions (the 75 in the Eg-type currents 
is harmless since it is attached to the heavy quark line). We shall not commit ourself 
to any definite 75-scheme but give our results for the two generic 75-schemes that are 
being used in the literature, namely the naive 75-scheme with an anticommuting 75 
(for a general exposition of the naive 75-scheme see ||23|| ) and the 't Hooft-Veltman- 
Breitenlohner-Maison (HV) 75-scheme In the naive 75-scheme one can simply 



anticommute the 75 in the AQ-currents to the end of the light fermionic string and 
then apply the general result Eq. (|33|). For the total anomalous dimension 7^ up to 
two loop order one then obtains (see Table 1 for the specific value of the pair (n, s)) 
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particles 


75 A ° 





+ 1 


Ai 


7 5 7o 


1 


-1 


A 2 


7 


1 


+1 


Si, 


7o7 


2 


-1 


S 2 , S2 


7 5 HV 


4 


-1 


Ai 


7 5 HV 7o 


3 


+1 


A 2 



Table 1: Specific values of the pair (n, s) for particular cases of the light-side Dirac 
structure 17 

^ = - 8 O + 5 (16C(2) + 40 *- 796) (s) 2 ' 

7A, = "4 (^)+i(16C(2) + 20^-322) (^-f , 

7E1 = - 4 (s) + 5 (16C(2) + 20 ^- 290) (s) 2 ' 
^-5 + > c(2) + 8 * + 324) 2 

The anomalous dimensions of the spin-3/2 currents Js*i and Js*2 coincide with 
those of the spin-1/2 currents Jsi and Js 2 , resp., because the respective light-side 
Dirac structures are identical. 

In the HV 75-scheme the anomalous dimensions of Jsi and Js 2 remain unchanged 
as remarked on before. Because one has 

75 = 7o7i7273 = ^i^a/^VVV 4 (48) 

in the HV 75-scheme, this particular definition of 75 is naturally included in the 
general representation of the 7-matrix string T used before in Eq. ( pi]) and Eq. (f43"|). 
In the HV 75-scheme one obtains 

7m = -8 (g) + i(16C(2) - 24iV F + 260) {^) , (49) 

7A2 = -4 (g) + 1(16C(2) - 12V F + 206) (^) 2 . (50) 

Comparing Eqs. fl44D and ( flBl) with Eqs. (||]) and (|50|) one observes that the ((2) 
contributions are not affected by the choice of 75, but there are big differences in the 
remaining contributions between the two 75-schemes. Needless to say that physical 
matrix elements do not depend on the choice of 75. 
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(44) 
(45) 
(46) 
(47) 



With the availability of the full two-loop result one can check on the "quick and 
easy" -prescription advanced in to estimate the magnitude of two-loop correc- 
tions. This prescription goes by the name of "naive non-abelianization" . The idea 
is to replace N F by N F — 33/2 in the easily computed light-quark loop contribution. 
This prescription is based on the observation that the first term of the /3-function is 
a large number and hence the dominant part of any corrections is connected with 
the term proportional to (3^ . Thus one can hope to arrive at a good estimate of 
a non-abelian result by replacing the abelian /J^-function by its non-abelian coun- 
terpart. We find that this prescription gives a good estimate of the full result with 
an accuracy of 15-30% as can easily be checked. An exception is the case Js 2 (and 
thereby Js*2) where the "naive non-abelianization" estimate is completely off the 
mark. But note that in the latter case the anomalous dimension is an order of 
magnitude smaller than in the other cases. 

As mentioned before, the anomalous dimensions of the baryonic currents J HV 
and J AC are different at two loop order. However, the currents can be related to 
each other by finite renormalization O, [Tjl , 



J AC (/i) = Z r J HV (/i). (51) 

The finite coefficients depend on the Dirac structure T, which occurs in the 
currents J = qTq and J = (q T CTq)Q in the mesonic case and baryonic case, re- 
spectively. The coefficients Zr may be determined by comparing the renormalized 
matrix elements in Eq. (|5T| ) or by calculating the logarithmic derivative of Eq. fl5"T|), 



_ = 7jHV - 7j ac = O r . (52) 

Let us introduce a factor A r through the relation Z r = 1 + A r (a s (//)/47r), where 
the coefficient can be calculated from 

Arf3(a s ) = -O r , (53) 

where /3(a g ) was defined in Eq. @. Using the results for the two-loop anomalous 
dimensions derived before one obtains 

16 in \ 2 

75 = -C B {l\C A -AN F T F )[^ , 

8 (a \ 2 

O 1Bli = O 15l0 = -C B (llC A -AN F T F )^j , (54) 

For completeness we have also included the Dirac structures T = 757^, 757*77- and 
757«7o in our discussion even though they do not correspond to the s-wave baryon 
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states discussed in this paper. By calculating Ay according to Eq. 
obtains 



one finally 



y 757i 



■^75 



y 7570 



B 



Air 

K'nl ^ 

47T 



(55) 



1. 



■^757i7j ^757i70 

Note that the finite renormalization coefficients for the baryonic currents differ from 



those in the mesonic case discussed in [13, WM by the replacement Cb — > Cp 



5 Conclusion 

We have calculated the anomalous dimensions of static heavy baryon currents at 
two-loop order. The anomalous dimensions are basic ingredients in the analysis of 
QCD sum rules in the heavy baryon sector. In order to improve on the accuracy of 
the analysis of the two-point and three-point sum rules in the literature one needs 
to avail of the next-to-leading order (NLO) radiative corrections to these sum rules. 
Part of the NLO corrections are determined by the two-loop anomalous dimensions 
of the baryon currents which have been calculated in this paper. We hope to return 
to the subject of the NLO baryonic sum rule corrections in the near future. 
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Appendix A 



In this appendix we list the coefficient functions of the heavy-light bare proper vertex 
Vq HV> up to two-loop order. We use the abbreviations E x = T(l — e)T(l + 2e) and 
E 2 = T(l — e) 2 r(l + At), where T(x) is Euler's Gamma-function. 



Ahi 
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h (hi 

°3 
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h (hl 
°9 



■ (hi 



b ( h i 



-2aE x 
D-A 

AE 2 



2 (3D 2 -2AD + AA)E 2 



{D-Af (D — Q)(D — A) 3 (D — 3) 

AEf A(D-5)(D-2)E 2 



(D-A) 3 (D-3) (D - 6)(D - Af{D - 3) 
AEo 



(D-A) 2 (D-3) 

2{D - 2)E 2 
D - 6)(D - A) 2 (D - 3) 

-8E 2 AE 2 



(D-A) 2 (D-3) (D — A) 2 (D — 3) 

-2{D - 2)E 2 
(D — 6)(D — A) 2 (D — 3) 







Eo 



(D -6)(D -A) 2 {D -3) 
—E 2 

(D -6)(D -A) 2 (D -3) 
-AE 2 



2DEo 



li — 



(D-A) 2 (D-3) (D - 6)(£> - A) 2 (D - 3) 

-A{D - l)E 2 
(D -G)(D -A) 2 (D -3) 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 

(A7) 

(A8) 
(A9) 

(A10) 

(All) 

(A12) 
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Appendix B 



In this appendix we list the coefficient functions of the light-light bare proper vertex 
Vq IV> up to two-loop order using the abbreviations Qi = T(l — e) 2 r(l + e)/T(l — 2e) 
and Q 2 = T(l - e) 3 r(l + 2e)/T(l - 3e). F(x) is Euler's Gamma- function. 

&0 '° " (i a) (D-4)(D-3) 

Ail) _ -Qi (m) 
^ " 2(£>-3) ( j 

/,(») zQi 

°' 2 ~" 2(D-4)(D-3) 

,(«) 4QD - 2)Q 2 4(/J-2)(/J 3 -9D 2 + 20/J-4)Q 2 



(£> - 4) 3 (/J - 3)(D - 1) (D - 6)(D - 4) 3 (/J - 3)(£> - 1)(3D - 10) 
, ( m (-D 3 - 12D 2 + 52D - 72)Q 2 (D 4 - 13D 3 + 6QD 2 - 1Q8D + 176)Q 2 

°1,2 = 



2(D - Af(D - 3) 2 (D - 1) (D — Q)(D - A) 3 (D - 3)(D - 1)(3D - 10) 

»,(") = 91 ?Ql (no) 

1.3 (D-4)( J D-3) 2 (,D- 1) (D-4)(D-l)(3D-10)(3D-8) K ' 

All) Ql Q2 

1.4 ~ 2(D-4) 2 (D-3) 2 (D-1) (D-4) 2 (D-l)(3D-10)(3D-8) 



All) 

°2,0 



2(D-2)(D 2 - 12D + 24)Q 2 
(D -A) 3 (D -3) 2 (D - 1) 

8(D - 2)(Z> 3 - 8£> 2 + 12D + 8)Q 2 
~ (D - 6%D - A) 3 (D -3)(D- I) (3D - 10) 
ah) 4(3D 2 - 18g + 26)Q 2 

2,1 (/J - 4) 3 (/J - 3) 2 (D - 1) 

4(5L> 5 - 71D 4 + 350L> 3 - 680£> 2 + 296D + 320)Q 2 
+ (D - Q)(D - 4) 3 (D -3)(D- 1)(3D - 10) (3D - 8) 
(£> 2 - 18D + 40)Q 2 



6 (H) _ 
"2,2 — 



(D — 4) 3 (D — 3) 2 (D — 1) 

4(4L> 4 - 43L> 3 + 12AD 2 - AOD - 160)Q 2 
+ (D - Q)(D - A) 3 (D -3)(D- 1)(3D - 10) (3D - 8) 

6 <«> = ^2 fB3) 

2,3 (D-4) 2 (D-3) 2 (D-1) (D-4) 2 (D-l)(3D-10)(3D-8) { ' 

All) Ql 2DQ 2 
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Appendix C 

In this appendix we present the two-loop coefficient functions of the irreducible 
heavy-light-light proper vertex Vq T \ 
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Figure Captions 



Fig. 1: One-loop (diagram i =0) and two-loop (diagrams i =1-11) heavy-light 
vertex corrections. Light lines: light quarks; heavy lines: static heavy 
quark; curly lines: gluons; dashed lines: ghost contributions 

Fig. 2: Two- loop heavy-light-light irreducible vertex correction. Light lines: 
light quarks; heavy lines: static heavy quark; curly lines: gluons. The 
remaining four diagrams are obtained by reflection on the heavy quark 
line 
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